DISTRIBUTION  STATEMENT  A:  Approved  for  public  release;  distribution  is  unlimited. 


High  Wind  Upper  Ocean  Mixing  with 
Explicit  Surface  Wave  Processes 

Peter  P.  Sullivan 

National  Center  for  Atmospheric  Research 
Boulder,  CO  80307-3000 

Phone:  (303)  497-8953  fax:  (303)  497-8171  email:  pps@ucar.edu 

James  C.  McWilliams 
Department  of  Atmospheric  Sciences  and 
Institute  of  Geophysics  and  Planetary  Physics,  UCLA 
Los  Angeles,  CA  90095-1565 

Phone: (310)  206-2829  fax: (310)  206-5219  email:  jcm@atmos.ucla.edu 
Grant  Number:  N00014-09-C-0088 


LONG-TERM  GOALS 

The  work  described  here  supports  the  Office  of  Naval  Research  Departmental  Research  Initiative 
(DRI)  for  research  on  the  “Impact  of  Typhoons  on  the  Ocean  in  the  Pacific”  (ITOP).  One  of 
central  themes  of  the  DRI  is  to  better  characterize  and  predict  the  ocean  boundary  layer  (OBL) 
and  its  impact  on  typhoon  (hurricane)  evolution.  This  is  one  component  of  developing  improved 
prediction  models  for  the  coupled  atmosphere-ocean-wave  system.  Cooling  of  the  sea  surface 
temperature  (SST)  is  a  critical  coupling  variable  influencing  atmosphere-ocean  hurricane  dy¬ 
namics;  SST  is  largely  determined  by  OBL  turbulence,  surface  wave  processes,  and  mixed  layer 
entrainment.  Our  research  goal  is  to  model  the  strongly  forced  wind  and  wave  driven  upper 
OBL  using  turbulence  resolving  large-eddy  simulation  (LES)  with  explicit  wave  effects,  viz., 
wave-current  interactions  and  breaking  waves  and  examine  their  impact  on  ocean  mixing  dur¬ 
ing  hurricane  events. 


OBJECTIVES 

The  specific  research  objectives  for  ITOP  are:  (1)  conduct  process  studies  using  LES  of  the 
OBL  with  different  combinations  of  time  varying  large  scale  forcings  and  surface  wave  effects 
and  examining  their  impact  on  ocean  mixing;  (2)  evaluate  and  compare  these  LES  results  with 
predictions  obtained  using  a  1-D  column  model  of  the  OBL  based  on  the  K-Prohle  Parameteri¬ 
zation  (KPP);  and  (3)  compare  our  simulation  results  with  available  observations.  Inertial  res¬ 
onance,  storm  residence  time,  and  the  larger  scale  environment  are  some  of  the  processes  to  be 
examined  in  our  simulations. 
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APPROACH 


In  order  to  examine  the  dynamics  of  the  upper  ocean  boundary  layer  under  high  wind  forcing 
we  developed  an  LES  model  of  the  OBL  that  accounts  for  surface  wave  effects  through  non¬ 
conservative  breaking  waves  and  phase-averaged,  conservative  wave-current  interactions  which 
lead  to  Langmuir  circulations  (Sullivan  et  al.  2007).  The  equations  for  the  resolved  flow  compo¬ 
nents  are  the  Craik-Leibovich  (CL)  theory  with  crucial  wave-current  coupling  through  the  vor¬ 
tex  force.  The  larger  scale  forcing  (momentum  and  scalar  fluxes)  and  wave  fields  are  externally 
imposed  in  the  LES. 


WORK  COMPLETED 

Meetings:  We  attended  a  PI  meeting  in  Santa  Fe,  NM  that  focused  on  a  first  look  at  the  ex¬ 
tensive  observation  datasets  collected  during  the  ITOP  held  campaign  conducted  in  the  Fall  of 
2010. 

OBL  LES  and  Coupling  to  Wavewatch  3:  During  the  past  year,  we  collaborated  exten¬ 
sively  with  Dr.  Leonel  Romero  from  Scripps  Institute  of  Oceanography,  who  helped  us  initiate  a 
one-way  coupling  between  the  wave  prediction  model  NOAA/NCEP  Wavewatch  3  and  our  LES 
code  for  the  OBL.  This  coupling,  although  very  idealized,  allows  us  to  drive  ocean  LES  with 
realistic  time  and  space  varying  wave  fields  that  are  generated  by  hurricane  winds.  Thus  a  con¬ 
sistent  set  of  wind  and  wave  fields  are  input  to  the  LES.  Our  canonical  surface  winds  and  waves 
are  an  idealization  of  Hurricane  Frances  (Sanford  et  al,  2011;  Zedler  2007).  In  addition  to  the 
standard  suite  of  variables  output  by  Wavewatch  3  (WW3),  the  wave  spectra  are  further  pro¬ 
cessed  to  generate  fields  pertinent  to  our  air-sea  interaction  research,  viz.,  wave  age,  turbulent 
Langmuir  number,  Stokes  drift  velocity,  wind-wave  alignment,  and  Stokes  depth  scale.  Prelimi¬ 
nary  results  from  the  WW3  simulations  are  discussed  in  the  RESULTS  section  and  will  be  fur¬ 
ther  described  in  future  publications. 

The  general  design  of  the  LES  experiments  is  to  pass  a  hurricane  vortex  over  small  LES  do¬ 
mains  as  sketched  in  Fig.  1.  Now,  the  surface  momentum  and  scalar  fluxes  and  wave  fields  im¬ 
posed  at  the  water  surface  are  those  generated  by  WW3  along  selected  time  slices  within  the 
storm  path.  Our  process  studies  examine  the  impact  of  inertially  rotating  wind  stress,  resonant 
and  anti-resonant  forcing,  along  lines  to  the  right  and  left  of  the  storm  track. 

Two  LES  computational  domains  with  different  horizontal  extent  are  used.  On  the  left  hand 
side  X  =  —60  km,  which  has  slightly  weaker  forcing,  a  small  domain  (Lx,  Ly,  Lz)  =  (750,  750,  —240) 
m  was  found  acceptable.  On  the  right  hand  side  of  the  storm,  X  =  60  km,  the  vigorous  forcing 
and  inertially  rotating  wind  stress  generates  a  system  of  large-scale  random  gravity  waves  in  the 
hurricane  wake.  Test  runs  indicated  that  these  structures  were  only  adequately  captured  by  in¬ 
creasing  the  horizontal  area  of  the  small  domain  LES  by  a  factor  of  4.  Thus,  a  large  domain  is 
used  on  the  right  hand  side  with  (Lx,  Ly ,  Lz)  =  (1500, 1500,  —240)  nr.  The  number  of  gridpoints 
in  the  small  domain  is  (Nx,  Ny,  Nz)  =  (512,  512,  256)  and  is  increased  to  (1024, 1024,  256)  in  the 
large  domain;  thus  the  horizontal  spacing  Ax  =  Ay  =  1.46  m  in  all  the  computational  boxes.  In 


the  vertical  direction,  we  use  a  stretched  vertical  grid  where  the  spacing  between  any  neighbor¬ 
ing  cells  is  held  fixed  at  the  ratio  K  =  1.00451  with  the  first  gridpoint  located  at  A z\  =  —0.5 
m.  Based  on  these  choices,  the  cell  aspect  ratio  Ax/ A z  =  (2.93, 1.54)  at  z  —  (0,  —100)  m.  The 
small  and  large  computational  boxes  are  able  to  capture  small-scale  turbulence  at  early  times 
and  larger  scales  of  motion  that  develop  during  the  rapid  deepening  of  the  OBL  at  the  time  of 
maximum  winds. 

The  initial  temperature  sounding  is  simply  constructed  from  two  linear  segments:  Over  the 
depth  hi  <  z  <  0,  6  —  302.4  K  and  for  0  <  hi  the  temperature  decreases  at  a  rate  of  0.04 
K  m-1.  This  is  approximately  the  variation  observed  for  Hurricane  Frances  (Sanford  et  al. , 
2011).  hi  =  —32  m  is  the  initial  OBL  depth.  The  Coriolis  parameter  is  /  =  0.681  x  10~4  s^1 
which  corresponds  to  an  inertial  period  T  =  25.6  hrs.  The  simulations  are  initiated  from  rest 
with  small  random  perturbations  in  the  temperature  field  under  the  wind  and  wave  forcing  from 
hour  10  of  the  WW3  simulation.  At  this  time,  the  wind  and  wave  fields  are  slowly  varying  with 
winds  less  than  10  ms-1.  Next,  these  external  forcings  are  held  constant  and  the  simulations  are 
integrated  for  10  physical  hours  to  generate  fully  developed  equilibrium  turbulence.  The  result¬ 
ing  mean  currents  are  0(0.05)  ms-1,  and  the  initial  state  of  these  weak  currents  has  minimal 
impact  on  the  resulting  long  time  ocean  response  under  strong  forcing  (Crawford  et  al,  1996, 
p.  889).  The  time  step  At  in  the  code  is  dynamically  picked  to  obey  a  fixed  Courant-Fredrichs- 
Lewy  condition.  At  varies  by  an  order  of  magnitude,  from  [0.3  —  5]  s,  and  more  than  300,000 
timesteps  are  needed  to  cover  the  70  physical  hours  of  the  simulation.  With  a  storm  translation 
speed  of  5.5  m  s”1  the  LES  domains  are  initially  located  more  than  800  km  from  the  storm  cen¬ 
ter  shown  in  Fig.  1.  The  parallelization  of  the  code  is  described  by  Sullivan  &  Patton  (2011). 

RESULTS 

The  X  —  t  Hovmollcr  contour  maps  presented  in  Figs.  2  and  3  illustrate  some  of  the  variability 
in  the  wind  and  wave  fields  induced  by  our  hypothetical  storm  over  the  horizontal  domain  X  = 
±300  km  for  a  time  period  of  more  than  50  hrs.  In  this  figure,  time  is  referenced  to  the  time  of 
maximum  winds  tm.  Thus,  the  X  cross  section  is  ahead  of  the  storm  for  t  —  tm  <  0  and  is  behind 
the  storm  (in  the  storm  wake)  for  t  —  tm  >  0.  Inspection  of  the  contours  in  Fig.  2  shows  that  the 
windspeed  |Ua|  is  asymmetrical  about  X  —  0  km,  a  consequence  of  the  storm  translation  speed 
Vh,  with  the  largest  values  ~  50  ms-1  on  the  right  hand  side  of  the  storm  track. 

Given  the  surface  winds  the  primary  output  of  the  wave  model  is  the  2D  wave  height  spectrum 
F(u 7,  <L)  which  varies  with  radial  frequency  w  and  angular  position  <f>.  The  wave  spectra  are 
further  processed  to  generate  the  wave-age  field  (Fig.  2)  and  Stokes  drift  velocity  (Fig.  3).  Wave 
age  is  defined  as  A  =  Cp/\XJa\  where  Cp  is  the  phase  speed  of  the  peak  in  the  wave  height  spec¬ 
trum.  A  varies  widely  in  space  and  time  across  the  storm  path.  Wave  age  is  a  broad  indicator 
of  wave  breaking  activity;  the  breaking  activity  tends  to  increase  with  decreasing  wave  age,  but 
also  depends  on  wave  steepness.  Extensive  periods  of  growing  waves  A  ~  0.5  are  observed  near 
the  time  of  maximum  winds  with  the  lowest  wave  age  on  the  right  hand  side  of  the  storm.  For 
time  periods  well  ahead  of  the  storm  center,  A  ~  1.2  indicative  of  near  wind- wave  equilibrium 
conditions.  At  t  —  tm  ~  —20  hrs,  the  wave  age  jumps  abruptly,  A  >  2,  as  fast  moving  swell  prop¬ 
agates  rapidly  ahead  of  the  storm  eye.  In  the  wake  of  the  storm  t  —  tm  >  0,  the  wave-age  again 


slowly  approaches  wind-wave  equilibrium  and  at  long  times  t  —  tm  >  20  hrs  the  wave  state  again 
becomes  swell  dominated.  These  estimates  of  wave  age  in  a  hurricane  are  consistent  with  val¬ 
ues  deduced  from  observations  of  ID  spectra  documented  by  Young  (1998)  who  reports  values 
ranging  from  [0.5  -  2.0]. 

A  key  variable  for  our  LES  is  the  Stokes  drift  velocity  us  which  appears  in  the  vortex  force  us  x 
uj  where  U)  is  the  resolved  vorticity.  For  a  broadband  wave  spectrum  us  is  obtained  by  the  inte¬ 
gration  (Kenyon,  1969) 

POO  P7V 

us(z)  =  2  /  /  4>)  e2^z  d(j)dw  (1) 

Jo  J-n 

where  k  denotes  the  horizontal  wavenumber  vector.  The  Stokes  drift  velocity  us  inherits  the 
complexity  of  the  wave  held  through  (1);  it  varies  with  wind  speed,  spatial  position,  and  time  as 
depicted  in  Fig.  3.  For  our  particular  choice  of  storm  propagation  direction,  the  largest  Stokes 
drift  occurs  in  component  vs  j  on  the  right  hand  side  of  the  storm  track.  The  distribution  of  the 
wave  held  is  more  asymmetrical  than  the  wind  held  mainly  due  to  the  “extended”  fetch  on  the 
right  hand  side  of  a  translating  hurricane  (Young,  2003).  As  a  result,  the  distribution  of  vs  is 
strongly  asymmetrical  about  X  —  0  km  a  reflection  of  the  strongest  winds  and  better  wind- 
wave  alignment  on  the  storm  right  hand  side.  us  i  is  clearly  largest  (most  negative)  in  front  of 
the  storm  center  with  contours  roughly  symmetrical  about  X  —  0  km. 

Past  LES  of  steady  aligned  wind-wave  cases  and  a  limited  number  of  held  observations  find  ver¬ 
tical  velocity  is  a  robust  signature  of  Langmuir  turbulence.  In  the  Langmuir  turbulence  regime, 
coherent  Langmuir  cells  form  and  these  structures  can  be  detected  in  the  instantaneous  vertical 
velocity  held  and  furthermore  impact  the  statistical  (average)  moments  of  w.  Figure  4  compares 
the  maximum  resolved  second  order  moment  ( w2)max  and  the  minimum  resolved  third-order  mo¬ 
ment  (w3)min  diagnosed  from  our  simulations  with  and  without  wave  effects  along  the  tracks 
X  =  ±60  km.  These  temporally  varying  moments  are  obtained  by  hrst  computing  the  vertical 
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prohles  ( w(z,t )  )  and  ( w(z,t )  )  at  all  z  for  a  given  t:  (  )  denotes  an  x  —  y  spatial  average  over 
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the  LES  domain.  (w(t)  )max  and  (w(t)  )min  are  found  by  searching  for  the  appropriate  maxi¬ 
mum  or  minimum  along  their  respective  vertical  prohles  and  are  then  normalized  by  the  local 
surface  friction  velocity  w*(A",  t).  We  note  the  vertical  locations  of  the  maximum  and  minimum 
vary  with  time.  These  positions  are  always  located  in  the  upper  surface  layer  of  the  OBL  but 
are  found  deeper  in  the  water  as  the  winds  and  waves  increase. 

The  time  series  shown  in  Fig.  4  are  compelling  evidence  for  Langmuir  turbulence  in  the  hurri¬ 
cane  driven  OBL.  The  w— moments  are  sensitive  to  the  location  within  the  storm  track,  depend 
on  the  time  and  space  evolving  wind  and  wave  helds,  and  the  inclusion  of  phase  averaged  wave- 
current  interactions  in  the  LES.  Over  the  lifetime  of  the  simulation  with  vortex  force  (th2)max  > 

2  along  the  track  X  =  60  km.  Notice  also  that  during  the  time  period  [-5,5]  hours,  the  max¬ 
imum  vertical  velocity  variance  exceeds  2.5  and  this  level  persists  in  the  hurricane  wake,  w* 
varies  by  about  a  factor  of  7  over  the  simulation  lifetime  and  thus  the  dynamic  range  of  (w2)max 
is  broad.  Surprisingly,  simple  scaling  by  w*  captures  much  of  the  variation  in  w— variance.  On 
the  left  side  of  the  storm  track,  the  maximum  w— variance  at  hrst  tracks  its  counterpart  on 
the  right  hand  side,  then  gradually  decays,  reaches  a  minimum  value  at  the  time  of  maximum 
winds,  and  gradually  recovers  in  the  hurricane  wake.  The  differences  between  the  variations  on 


the  right  and  left  hand  side  of  the  storm  are  dramatic  and  can  be  traced  to  the  wave  fields  and 
wind-Stokes  alignment.  The  mis-aligned  wind  and  wave  fields  and  the  very  low  values  of  Stokes 
drift  around  the  time  of  maximum  winds  on  the  left  hand  side  of  the  storm  track  effectively  shut 
off  the  vortex  force.  In  the  simulations  without  vortex  force,  the  maximum  w— variance  is  nearly 
independent  of  the  track  position  and  is  0(1),  similar  to  what  is  expected  for  a  flat-wall  bound¬ 
ary  layer  flow. 

The  impact  of  wave-current  interactions  on  the  third-order  vertical  velocity  moment  is  even 
more  impressive  upon  comparing  cases  with  and  without  vortex  force.  With  vortex  force  the 
minimum  third-order  moment  can  be  10  times  larger  than  its  no  vortex  force  counterpart.  Thus 
we  find  CL  mechanics  are  active  and  critical  in  our  simulations  but  varies  with  the  wave  state. 

The  surface  Stokes  drift  velocity  on  the  right  side  of  the  storm  track  scales  reasonably  well  with 
the  local  friction  velocity  w*.  Our  results  are  in  good  agreement  with  the  measurements  of  D’Asaro 
(2001). 

Flow  visualization  of  the  instantaneous  vertical  velocity  at  the  time  of  maximum  winds  is  com¬ 
pelling  evidence  for  Langmuir  turbulence.  In  Fig.  5,  we  show  images  of  resolved  w/u *  from  two 
simulations:  the  upper  panel  is  from  the  left  (non-resonant)  side  while  the  lower  panel  is  from 
the  right  (resonant)  side  of  the  storm  track.  As  expected  based  on  the  time  series  of  (w2)max 
and  ( w3)min  shown  in  Fig.  4,  there  is  scant  evidence  for  the  formation  of  Langmuir  cells  on  the 
left  hand  side  of  the  track.  In  contrast,  impressive  vertically  oriented  streaks  are  found  on  the 
right  hand  side.  Note  at  this  time  the  winds  and  stresses  in  the  lower  image  are  oriented  in  the 
positive  y  direction  (see  vectors  in  the  image).  The  streaks  show  mergers  and  are  found  to  sig¬ 
nificantly  alter  the  turbulent  dynamics  in  the  upper  OBL.  Thus,  wave-current  interactions  are 
most  active  on  the  resonant  side  of  the  storm. 


IMPACT  /APPLICATIONS 

The  LES  results  obtained  here  for  hurricane  driven  OBLs  can  be  used  to  guide  the  interpreta¬ 
tion  of  observations  collected  during  the  ITOP  program.  In  addition,  the  results  can  be  used  to 
test  simpler  1-D  parameterizations  of  the  ocean  mixed  layer  that  are  used  in  large  scale  models. 
Our  particular  interest  is  in  evaluating  and  improving  the  so-called  K-prohle  parameterization 
(KPP)  that  is  routinely  used  in  the  Regional  Ocean  Modeling  System  (ROMS). 


TRANSITIONS  &  RELATED  PROJECTS 

The  present  work  has  links  to  the  ONR  DRI  on  High  Resolution  Air-Sea  Interaction  (HIRES) 
that  focuses  on  the  interaction  of  waves  and  turbulence  in  the  atmospheric  surface  layer.  For 
a  description  of  HRES  see  http://airsea.ucsd.edu/hires/.  The  LES  model  being  developed  for 
HIRES  is  also  being  used  in  the  present  work. 


REFERENCES 


D’Asaro,  E.,  2001:  Turbulent  vertical  kinetic  energy  in  the  ocean  mixed  layer.  Journal  of  Physi¬ 
cal  Ocea?iography,  31,  3530-3537. 

Grant,  A.  L.  M.  &  S.  E.  Belcher,  2009:  Characteristics  of  Langmuir  turbulence  in  the  ocean 
mixed  layer.  Journal  of  Physical  Oceanography ,  39,  1871-1887. 

Kenyon,  K.  E.,  1969:  Stokes  drift  for  random  gravity  waves.  Journal  of  Geophysical  Research , 

74,  6991-6994. 

McWilliams,  J.C.,  P.  P.  Sullivan  &  C-H.  Moeng,  1997:  Langmuir  turbulence  in  the  ocean.  Jour¬ 
nal  of  Fluid  Mechanics, 334,  1-30. 

Sanford,  T.  B.,  J.  F.  Price  &  J.  B.  Girton,  2007:  Highly  resolved  observations  and  simulations  of 
the  ocean  response  to  a  hurricane.  Journal  of  Physical  Oceanography,  41,  1041-1056. 

Sullivan,  P.P.,  J.C.  McWilliams,  &  W.K.  Melville,  2007:  Surface  gravity  wave  effects  in  the 
oceanic  boundary  layer:  Large-eddy  simulation  with  vortex  force  and  stochastic  breakers.  Jour¬ 
nal  of  Fluid  Mechanics,  593,  405-452. 

Sullivan,  P.  P.  &  E.  G.  Patton,  2011:  The  effect  of  mesh  resolution  on  convective  boundary  layer 
statistics  and  structures  generated  by  large-eddy  simulation.  Journal  of  the  Atmospheric  Sci¬ 
ences,  [in  press]. 

Tseng,  R-S.  &  E.  D’Asaro,  2004:  Measurements  of  turbulent  vertical  kinetic  energy  in  the  ocean 
mixed  layer  from  Lagrangian  floats,  Journal  of  Physical  Oceanography,  34,  1984-1990. 

Young,  I.  R.,  1998:  Observations  of  the  spectra  of  hurricane  generated  waves.  Ocean  Engineer¬ 
ing,  25,  number  4-5,  261-276. 

Young,  I.  R.,  2003:  A  review  of  the  sea  state  generated  by  hurricanes.  Marine  Structures,  16, 
201-218. 

Zedler,  S.E,  2007:  Strong  wind  forcing  of  the  ocean.  PhD  Thesis,  University  of  California,  San 
Diego,  San  Diego,  California,  pp.  127. 


PUBLICATIONS 

Sullivan,  P.  P.  &  E.  G.  Patton,  2011:  The  effect  of  mesh  resolution  on  convective  boundary  layer 
statistics  and  structures  generated  by  large-eddy  simulation.  Journal  of  the  Atmospheric  Sci¬ 
ences,  [in  press]. 

Weil,  J.  C.,  P.P.  Sullivan,  E.  G.  Patton  &  C-H.  Moeng,  2011:  Statistical  variability  of  dispersion 
in  the  convective  boundary  layer:  LPDM-LES  model  ensembles  and  observations.  Boundary- 
Layer  Meteorology,  [submitted]. 

Nilsson,  E.  O.,  A.  Rutgersson  &  P.  P.  Sullivan,  2011:  Convective  boundary  layer  structure  in 
the  presence  of  wind- following  swell.  Quarterly  Journal  of  the  Royal  Meteorological  Society, 
[submitted] . 


Van  Roekle,  L.  P.,  B.  Fox-Kemper,  P.  P.  Sullivan,  P.  E.  Hamlington  &  S.  R.  Haney,  2011:  The 
form  and  orientation  of  Langmuir  cells  for  misaligned  winds  and  waves.  Journal  of  Geophysical 
Research ,  [submitted]. 

Hanley,  K.  E.,  S.  E.  Belcher  &  P.  P.  Sullivan,  2011:  Response  to  “Comments  on  a  global  clima¬ 
tology  of  wind- wave  interaction”.  Journal  of  Physical  Oceanography ,  [in  press]. 

Moeng,  C.-H.  and  P.P.  Sullivan  2011:  Large  Eddy  Simulation.  Encyclopedia  of  Atmospheric 
Sciences  2nd  Edition,  Eds.  G.  North,  F.  Zhang  and  J.  Pyle.  Academic  Press,  [submitted]. 

Liang,  J.,  J.  C.  McWilliams,  P.  P  .  Sullivan  &  B.  Baschek,  2011:  Modeling  bubbles  and  dis¬ 
solved  gases  in  the  ocean.  Journal  of  Geophysical  Research  -  Oceans.  116,  C03015,  1-17. 

Kukulka,  T.,  A.  J.  Plueddemann,  J.  H.  Trowbridge,  &  P.  P.  Sullivan,  2011:  The  influence  of 
crosswind  tidal  currents  on  Langmuir  circulation  in  a  shallow  ocean.  Journal  of  Geophysical  Re¬ 
search  - Oceans ,  116,  C08005,  1-15. 

Lenschow,  D.  H.,  M.  Lothon,  S.  D.  Mayor,  P.  P.  Sullivan,  &  G.  Canut,  2011:  A  comparison  of 
higher-order  vertical  velocity  moments  in  the  convective  boundary  layer  from  lidar  with  in  situ 
measurements  and  LES.  Boundary-Layer  Meteorology ,  [in  press]. 

Suzuki,  N.,  T.  Hara,  &  P.  P.  Sullivan,  2011:  Turbulent  airflow  at  young  sea  states  with  frequent 
wave  breaking  events:  Large  eddy  simulation.  Journal  of  the  Atmospheric  Sciences ,  68.  1290- 
1305. 

Patton,  E.  T.  Horst,  P.  Sullivan,  D.  Lenschow,  S.  Oncley,  W.  Brown,  S.  Burns,  A.  Guenther,  A. 
Held,  T.  Karl,  S.  Mayor,  L.  Rizzo,  S.  Spuler,  J.  Sun,  A.  Turnipseed,  E.  Allwine,  S.  Edburg,  B. 
Lamb,  R.  Avissar,  R.  Calhoun,  J.  Kleissl,  W.  Massman,  K.  Paw-U,  &  J.  Weil,  2011:  The  canopy 
horizontal  array  turbulence  study  (CHATS).  Bulletin  of  the  American  Meteorological  Society , 

92,  593-611. 

Moeng,  C.-H.  Moeng,  P.  P.  Sullivan,  M.  F.  Khairoutdinov,  &  D.  A.  Randall,  2010:  A  mixed 
scheme  for  subgrid-scale  fluxes  in  cloud-resolving  models.  Journal  of  the  Atmospheric  Sciences , 

67,  3692-3705. 

Sullivan,  P.  P.  and  E.  G.  Patton,  2011:  Offshore  marine  boundary-layer  winds  predicted  by  a 
large  eddy  simulation  model  with  resolved  surface  waves.  6th  Theoretical  Fluid  Mechanics  Con¬ 
ference,  American  Institute  of  Aeronautics  and  Astronautics,  Honolulu,  HI. 

Nguyen,  K.  X.  S.  P.  Oncley,  T.  W.  Horst,  P.  P.  Sullivan  &  C.  Tong,  2010:  Investigation  of  subgrid- 
scale  turbulence  in  the  atmospheric  surface  layer  using  AHATS  held  data.  American  Physical 
Society,  Fluid  Dynamics ,  Long  Beach,  CA. 

Suzuki,  N.,  T.  Hara,  &  P.  P.  Sullivan,  2010:  Turbulent  airflow  at  young  sea  states  with  frequent 
wave  breaking  events:  Large  eddy  simulation.  17th  Conference  on  Air  Sea  Interaction ,  Annapo¬ 
lis,  MD. 

Liang,  J.,  J.  C.  McWilliams,  P.  P  .  Sullivan,  2010:  Modeling  gas  bubbles  in  the  ocean  boundary 
layer.  17th  Conference  on  Air  Sea  Interaction,  Annapolis,  MD. 


t 


t  +  5t 


Vh 


Figure  1:  Sketch  of  an  idealized  hurricane  vortex  used  for  setting  the  surface  forcing  conditions  in 
LES.  The  vortex  propagates  left  to  right  with  speed  Vh  and  has  a  characteristic  length  scale  (radius 
of  maximum  winds)  A.  The  families  of  solid  and  dotted  circles  indicate  the  position  of  the  vortex  at 
initial  time  t  and  at  a  later  time  t  +  At,  respectively.  The  fixed  points  (R,  L)  to  the  right  ( X  >  0)  and 
left  (X  <  0)  of  the  vortex  center  (c)  feel  the  time  history  of  winds  and  waves  along  the  horizontal  gray 
lines. 
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Figure  2:  The  hurricane  vortex  propagates  left  to  right  and  the  fields  are  shown  in  the  form 
of  an  X-t  Hovmoller  diagram,  i.e.,  the  temporal  variation  of  a  particular  held  is  shown  for  a 
chosen  x  cross  section.  Time  is  referenced  to  the  time  of  maximum  winds  tm.  The  upper  panel 
shows  contours  of  wind  speed  (in  m  s-1J  imposed  in  the  Wavewatch  3  simulations.  The  lower 
panel  depicts  the  wave  age  A  =  C'p/|Ua  over  the  lifetime  of  the  storm. 
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Figure  3:  Upper  and  lower  panels  show  the  Stokes  velocity  components  (us,vs)  at  the  water 
surface  for  the  imposed  wind  field.  The  color  bar  is  in  units  of  ms4. 
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Figure  4:  Temporal  variation  of  the  maximum  resolved  vertical  velocity  variance  (w2)max  and 
minimum  resolved  third-order  moment  (w3)min  for  cases  with  and  without  wave  effects  and  lo¬ 
cations  left  and  right  of  the  storm  center.  Simulations  with  wave  effects  [a)  X  =60  km,  b)  X  = 
-60  km]  and  simulations  with  no  wave  effects  [c)  X  =  60  km,  d)  X  =  -60  km].  The  normaliza¬ 
tion  is  by  u,(X,  t). 


Figure  5:  Flow  visualization  of  the  vertical  velocity  field  normalized  by  local  u *  on  the  left 
(upper  panel )  and  right  (lower  panel )  sides  of  the  storm  at  the  time  of  maximum  winds.  The 
arrows  indicate  the  direction  of  the  wind  (white)  and  Stokes  drift  (red).  Note  the  formation  of 
vigorous  Langmuir  cells  on  the  right  hand  side  of  the  storm  track. 


